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(57) Magnetic resonance is excited (50) in first and 
second species dipoles of a subject in a temporally con- 
stant magnetic field. The resonance is refocused (52) to 
generate a spin echo (54) centered at a time when the 
first and second species resonance signals are in- 
phase. Gradients echoes (64, 68) are generated, cen- 
tered at a time (2n+1)7i/5o before and after the spin 
echo, where 5co is a difference between the first and sec- 
ond species resonance frequencies. In this manner, the 
first and second species signals are 180° out-of-phase 
in the gradient echoes. The resonance is refocused (82) 
one or more times to generate additional spin and gra- 
dient echoes with different phase encodings (78). The 
sequence is repeated with yet more phase encodings, 
and magnetic resonance signals from the spin echo and 
the two gradient echoes are reconstructed (86) into a 
spin echo image (s 0 ) and a pair of gradient echo images 
(s +1 , s.-,). A phase map is generated (90) from the spin 
and gradient echo images. One of the gradient echo im- 
ages is corrected (116) with the phase map. The phase 
corrected gradient image is additively combined (118) 
with the spin echo image to generate a first species im- 
age (112) and is subtractively combined (120) to gener- 
ate a second species image (114). 
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Description 

The present invention relates to a magnetic reso- 
nance imaging method and apparatus. It finds particular 
application in conjunction with the separation of water 
and fat signals from human patients and will be de- 
scribed with particular reference thereto. However, it is 
to be appreciated that the present technique is also ap- 
plicable to imaging sequences in which components ol 
the imaged region have very close resonance frequen- 
cies. 

Heretofore, subjects have been positioned in a tem- 
porally constant magnetic field such that selected di- 
poles preferentially align with the magnetic field. Radio 
frequency signals have been applied to cause magnetic 
resonance of the preferentially aligned dipoles. The ra- 
dio frequency magnetic resonance signal from the res- 
onating dipoles has been read out for reconstruction into 
an image representation. 

To strengthen the magnetic resonance signal, the 
resonance can be excited with a 90° radio frequency 
pulse followed by a 180° refocusing pulse. The 180° re- 
focusing pulse causes the resonating spin system to re- 
focus as a spin echo. The time between the 180° refo- 
cusing pulse and the spin echo is the same as the time 
between the 90° excitation pulse and the 180° refocus- 
ing pulse. 

Magnetic resonance echoes can also be induced 
by other disturbances of the spin system, such as re- 
versing the polarity of a magnetic field gradient to induce 
a gradient echo. 

As illustrated in U.S. Patent No. 4,833,407 of Hol- 
land, Provost, DeMeester, and Denison, spin and gra- 
dient echoes have been induced interleaved in the same 
repetition of the magnetic resonance imaging se- 
quence. Briefly summarized, an RF excitation and refo- 
cusing pulse were applied to induce a spin echo. The 
polarity of a magnetic field gradient along a read axis 
was reversed one or more times to induce one or more 
gradient echoes. When used with fast spin echo (FSE) 
techniques, refocusing radio frequency pulses are ap- 
plied after each spin echo to induce yet another spin 
echo. 

As indicated above, water and fat have close reso- 
nance frequencies, but differ by about 220 Hertz in a 
magnetic field of 1 .5 Tesla. In the 1 .5 Tesla field, the 220 
Hertz higher frequency component gains a fulJ revolu- 
tion on the slower component every 4.46 msec. That is, 
the signals are in-phase with a 4.46 msec, periodicity. 
By inducing one echo at 2.23 msec, after refocusing and 
another echo at 4.46 msec, after refocusing, a pair of 
echoes can be induced and their signals read out. One 
of the echo signals has the water and fat in-phase and 
the other has the water and fat signals 180° out-of- 
phase. The two echoes can be induced in two different 
repetitions of the spin echo imaging sequence by shift- 
ing the radio frequency refocusing pulse by 2.23 msec. 
See, for example, Glover, et al., "Three-Point Dixon 



Technique For True Water/Fat Decomposition with Bo 
Inhomogeneity Correction", Magnetic Resonance in 
Medicine, Vol. 1 8, pp. 371 -383 (1 991 ). One of the prob- 
lems with this technique is that the system is typically 
5 not sufficiently linear that one can add the in-phase and 
out-of-phase signals to get a water only signal and sub- 
tract the in-phase and out-of-phase signals to obtain a 
fat only signal. The time evolution ol the fat and water 
signals is sufficiently non-linear that in-phase and out- 
of-phase components do not cancel completely. To cor- 
rect for this non-linearity, the Glover, et al. technique 
generates three echoes - one at a nominal echo time, 
one 2.23 msec, advanced from the nominal echo time, 
and one 2.23 msec, retarded from the nominal echo 
time. Through the use of post-processing, Glover seeks 
to remove the non-linearities such that components of 
2.23 msec, displaced echoes substantially cancel. An- 
other drawback of the Glover technique is that three rep- 
etitions of the imaging sequence are required to gener- 
ate magnetic resonance echoes with retarded, ad- 
vanced, and reference timings. 

In accordance with one aspect of the present inven- 
tion, a magnetic resonance imaging method is provided. 
A resonance excitation pulse and a resonance refocus- 
ing pulse are applied to generate a spin echo. Read gra- 
dients are reversed in order to induce first and second 
gradient echoes in association with each spin echo. The 
spin echo is induced when the resonance signals of a 
first species and a second species are substantially in- 
phase. The gradient echoes are generated at a time 
when the resonance signals ol the first and second spe- 
cies are substantially 180° out-of-phase. In this manner, 
an in-phase spin echo and out-of-phase gradient ech- 
oes are generated in each repetition of the sequence. 

In accordance with another aspect ol the present 
invention, a magnetic resonance imaging apparatus is 
provided. A primary magnet generates a temporally 
constant magnetic field. A gradient magnetic field coil 
generates magnetic field gradients across an imaging 
region. One or more radio frequency coils induces and 
refocuses magnetic resonance of dipoles ol first and 
second species within the imaging region and receives 
magnetic resonance signals from the resonating first 
and second species dipoles in the examination region. 
The first and second species dipoles have resonance 
frequencies which differ by a frequency ya>. A timing and 
control circuit controls the radio frequency and gradient 
coils (1 ) to excite magnetic resonance, (2) to refocus the 
magnetic resonance to induce a spin echo, and (3) to 
induce gradient echoes centered at times nn/yco before 
and after a center of the spin echo, where n is an odd 
integer. A reconstruction processor (1) reconstructs 
magnetic resonance signals from the spin-echo into a 
spin echo image, (2) reconstructs magnetic resonance 
signals from the first gradient echo into a first gradient 
echo image, and (3) reconstructs magnetic resonance 
signals from the second gradient echo into a second 
gradient echo image. 
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One advantage of the present invention is that the 
reference, advanced, and retarded resonance echo sig- 
nals are collected in a single acquisition. 

Another advantage ot the present invention is that 
data acquisition is accelerated. 5 

Another advantage of the present invention is that 
it can be used to differentiate among more than two spe- 
cies. 

The invention will now be further described, by way 
of example, with reference to the accompanying draw- 
ings in which:- 

FIGURE 1 is a diagrammatic illustration of a mag- 
netic resonance imaging device in accordance with 
the present invention; 

FIGURE 2 is a diagrammatic illustration of a mag- 
netic resonance imaging pulse sequence in accord- 
ance with the present invention; and, 
FIGURE 3 is an alternate embodiment of the mag- 
netic resonance imaging pulse sequence. 

With reference to FIGURE 1 , a main magnetic field 
control 10 controls superconducting or resistive mag- 
nets 12 such that a substantially uniform, temporally 
constant magnetic field is created along a z-axis through 
an examination region 14. A magnetic resonance echo 
means applies a series of radio frequency (RF) and 
magnetic field gradient pulses to invert or excite mag- 
netic spins, induce magnetic resonance, refocus mag- 
netic resonance, manipulate magnetic resonance, spa- 
tially and otherwise encode the magnetic resonance, to 
saturate spins, and the like to generate magnetic reso- 
nance imaging and spectroscopy sequences. More spe- 
cifically, gradient pulse amplifiers 20 apply current puls- 
es to selected ones or pairs of whole body gradient coils 
22 to create magnetic field gradients along x, y, and z- 
axes of the examination region 14. A digital radio fre- 
quency transmitter 24 transmits radio frequency pulses 
or pulse packets to a whole body RF coil 26 to transmit 
RF pulses into the examination region. A typical radio 
frequency pulse is composed of a packet of immediately 
contiguous pulse segments of short duration which tak- 
en together achieve a selected magnetic resonance ma- 
nipulation. The RF pulses are used to saturate, excite 
resonance, invert magnetization, refocus resonance, or 
manipulate resonance in selected portions of the exam- 
ination region. For whole body applications, the reso- 
nance signals are commonly picked up by the whole 
body RF coil 26. 

For generating images of limited regions of the sub- 
ject, local coils are commonly placed contiguous to the 
selected region. For example, an insertable head coil 
30 is inserted surrounding a selected brain region at the 
isocenler of the bore. The insertable head coil optionally 
includes local gradient coils 32 which receive current 
pulses from the gradient amplifiers 20 to create magnet- 
ic field gradients along x, y, and z-axes in the examina- 
tion region within the head coil. A local radio frequency 



4 

coil 34 is used to excite magnetic resonance and receive 
magnetic resonance signals emanating from the pa- 
tient's head. Alternatively, a receive-only local coil can 
be used in conjunction with body-coil transmission. An 
RF screen 36 blocks the RF signals from the RF head 
coil from inducing eddy currents in the gradient coils and 
the surrounding structures. 

With continuing reference to FIGURE 1 and further 
reference to FIGURE 2, a sequence control circuit 40 
controls the gradient pulse amplifiers, the digital trans- 
mitter, and a digital radio frequency receiver 38 to gen- 
erate the pulse sequence of FIGURE 2. A timing means 
42 controls an Excitation RF pulse generator 44 and a 
retocusing RF pulse generator 46. The resonance exci- 
tation pulse generator 44 controls the radio frequency 
transmitter to cause a resonance excitation RF pulse 50. 
The refocusing RF pulse generator 46 causes the radio 
frequency transmitter to generate a refocusing pulse 52 
at a time t after the excitation pulse. The excitation pulse 
followed by a refocusing pulse causes a spin echo 54 
to occur a time 2t after excitation. The times t are se- 
lected such that the fat and water resonance signals are 
substantially in-phase at the center of the spin echo 54. 
That is, the time 2x is selected to be an integer multiple 
of resonance oscillations at both the resonance frequen- 
cy of water and the resonance frequency of fat, i.e., 2nn/ 
8tD, where 6(o is the frequency difference between the 
water and fat protons. Preferably, the time % is also se- 
lected such that the fat and water are substantially in- 
phase when the refocusing pulse 52 is applied. 

The timing means 42 also controls a read gradient 
control 60. The read gradient control 60 causes a first 
read gradient 62. The application of the read gradient 
62 is timed such that a gradient echo 64 is induced at a 
time (2n+1)n/8co before the center of the spin echo 54. 
A center lobe 66 of the read gradient is applied centered 
on the spin echo 54. The read gradient control again 
reverses to induce a second gradient echo 68 at time 
(2n+1)n/8a) after the spin echo and concurrently with a 
read gradient lobe 70. In this manner, a spin echo is gen- 
erated surrounded on either side by a pair of gradient 
echoes which are 180° or n out-of-phase with the spin 
echo. Again, for fat and water in a 1.5 Tesla field, W 
5co=2.23 msec. 

The timing means also controls a slice select gra- 
dient control 72 which controls the gradient amplifiers 
20 to cause slice select gradient pulses 74 concurrently 
with the RF excitation and refocusing pulses. The slice 
select gradients limit the examined region to a single 
slice. Alternately, slab select pulses may be applied and 
phase encode pulses can be applied to encode the re- 
sultant spin and gradient echoes along the slice select 
axis. The timing means further controls a phase encode 
gradient controller 76 which causes the gradient ampli- 
fiers 20 to induce phase encode gradient pulses 78 be- 
fore each of the spin and gradient echoes. The phase 
encode gradient echoes are applied such that the spin 
echo and its two associated gradient echoes have the 
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same phase encoding. Preferably, a dephasing gradient 
pulse 80 can be applied along the phase encode axis to 
remove or zero the phase encoding after the last gradi- 
ent echo. 

In the preferred embodiment, the radio frequency 
refocusing pulse controller 46 causes subsequent refo- 
cusing pulses 82 at a time 2t after the first refocusing 
pulse and at 2x intervals thereafter. The read, slice, and 
phase encode gradient controllers cause like gradient 
pulses following each refocusing pulse to collect a sec- 
ond and subsequent sets of image data. 

With reference again to FIGURE 1, the magnetic 
resonance signals from the spin echo and the pair of 
gradient echoes are received by the whole body RF coil 
26 or the localized coil 34 and conveyed to the digital 
receiver 36. A sorter 84 sorts the signals from the refer- 
ence, retarded, and advanced echoes. A reconstruction 
processor 86, preferably three parallel processors, re- 
constructs a reference image s 0 , a first, retarded image 
s. 1( and a second, advanced image s +1 . 
The images are defined by: 



s t = (P w - p F ) Q 



i(-4^H) 



s o = (Pw + PF) e 



: (P w " Pf) 



(1a), 



(1b), 



(1c). 



where p w and p F denote the real water and fat compo- 
nents, respectively, ^represents a complex phase to the 
local field error (also known as background phase), and 
<j> rf is the systematic phase error introduced from the RF 
hardware chain. The s. v s 0 , and s +1 images are stored 
in image memories 88., 88 0 , 88 + , respectively 

A phase map or error generator 90 analyzes the im- 
ages s. v s 0 , and s +1 to determine a map of the phase 
error at each pixel of an imaged slice. For calculational 
simplicity, a transform circuit 92 generates complex 
phase signals S + and S. defined by: 



^+3,1 ^ (P W -P F ) 
2s 0 (P W +P F ) 



cos$ 



S = —p. = -. : i sin$ 

2s 0 (P W +P F ) 



(2a), 



(2b). 



The complex phase due to the local field error, i.e., the 
background phase, is given by: 

Arg(Re(S + )+ i Im(S_) ) = 4 ± 2rin - np (3), 



where n denotes an integer and p varies between 0 and 
1. A phase map processor 94 processes the output of 
the transform 92 to generate a phase map in accord- 
ance with Equation (3) for storage in a phase map mem- 
ory 94. Of course, the phase map has both 2n and IT 
discontinuities. An unwrapping algorithm or processor 
96 unwraps or adjusts the phase map to eliminate the 
n and 2n discontinuities to recover the background 
phase. In the preferred embodiment, the background 
phase is fit to a polynomial <i> rrt using a least squares 
processing routine 98. In the preferred embodiment in 
which the background phase is fit to a polynomial, it is 
defined by: 
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♦fit = ^0 + a i x + b iV + a 2 x * + b 2v 2 + c 2*y + ■ 



(4). 



The background phase e* fil for each pixel is stored in a 
phase correction or background phase memory 100. 

A corrected image generator 110 combines the 
phase correction and the uncorrected reconstructed im- 
ages to generate phase corrected water and fat images 
which are stored in a water image memory 112 and a 
fat image memory 11 4. In the preferred embodiment, the 
water and fat images are defined by: 



water image - |s 0 +s^e 



fat image - 



(5a), 



(5b). 



More specifically, a multiplier 116 multiplies one of the 
s., and s +1 images with the phase correction from the 
background phase memory 100. An image adder 118 
adds the complex phase corrected gradient echo image 
with the complex spin echo image to generate the water 
image (Equation (5a)). A subtraction circuit 1 20 subtrac- 
tively combines the complex phase corrected gradient 
echo with the complex spin echo image to generate the 
fat image (Equation (5b)). Optionally, a weighting adjust- 
ment 122 is provided for multiplying the phase adjusted 
complex gradient echo image by an adjustable weight- 
ing factor to adjust the relative weighting between the 
spin echo image and the gradient echo image. 

A video processor 1 24 converts selectable slices or 
other portions of the water or fat images or combinations 
of the two into appropriate format for display on a mon- 
itor 126 or other human readable display device. 

Of course, this technique can also be used to dis- 
tinguish between other than fat and water. With refer- 
ence to FIGURE 3, the read gradient can have five lobes 
for distinguishing among three species. The values of t 
are selected such that the spin echo 54 occurs when all 
three species are in-phase. A first pair of gradient ech- 
oes 64, 68 for a first species are generated at a time n/ 
§co' before and after the spin echo, where 5oo' is the f re- 
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quency difference between the resonance frequency of 
the second species and the first species. A second pair 
of gradient echoes 130, 132 are generated a time nil/ 
5od' before and after the spin echo, where 5cd' is the fre- 
quency difference between the first and third species 
and n is an integer, preferably 1 or an odd number. The 
above-described phase correction process is repeated 
to generate a phase correction between the first and 
second species and another phase correction between 
the First and third species. The sum of the complex gra- 
dient echo image that is adjacent to the spin echo and 
the complex spin echo image itself provides an image 
of the first species. The sum of the complex gradient 
echo images that are farthest from the spin echo and 
the complex spin echo image itself provide an image of 
the first species, for example, water. The subtraction of 
the same sets of images provide images of the second 
and third species, for example, fat and silicone, respec- 
tively. 

Claims 

1. A magnetic resonance imaging method in which a 
resonance excitation pulse and a resonance refo- 
cusing pulse are applied to generate a spin echo 
(54) and in which read gradients are reversed in or- 
der to induce first and second gradient echoes (64, 
6B) in association with each spin echo, character- 
ized by: 

inducing the spin echo (54) when resonance 
signals of a first and second species are sub- 
stantially in-phase; and, 
generating the gradient echoes (64, 68) at a 
time when the resonance signals of the first and 
second species are substantially 180° out-of- 
phase such that an in-phase spin echo and out- 
of-phase gradient echoes are generated in 
each repetition of the imaging sequence. 

2. A method as claimed in claim 1 , including: 

reconstructing a spin echo image from the spin 
echo signals from a plurality of repetitions of the 
imaging sequence; 

generating first and second gradient echo im- 
ages from the gradient echo signals of the plu- 
rality of repetitions of the imaging sequence; 
generating a phase correction map from the 
spin echo image and the first and second gra- 
dient echo images; 

correcting phase error in the first gradient echo 
image in accordance with the phase map; 
additively combining the spin echo image with 
the phase corrected first gradient echo image 
to generate a first species image (112); and 
subtractively combining the spin echo image 



with the phase corrected first gradient echo im- 
age to generate a second species image (114). 

3. A method as claimed in claim 2, further including: 

5 

generating the first gradient echo (64) a time Ul 
5cd before the spin echo, where 6© is a differ- 
ence in the resonance frequencies of the first 
and second species; and, 
10 generating the second gradient echo (68) at a 

time n/5© after the spin echo. 

4. A method as claimed in claim 3, further including: 

1$ generating the spin echo (54) at a time when 

resonance signals of the first species : second 
species, and a third species are all in-phase; 
generating a third gradient echo (1 30) at a time 
n/5co' after the spin echo, where 5co' is a drfter- 

20 ence in the resonance frequencies between the 

first and third species; and 
generating a fourth gradient echo (132) at a 
time n/8®' before the spin echo. 

25 5. A method as claimed in any of the preceding claims, 
wherein the first species (11 2) is water and the sec- 
ond species (114) is fat. 

6. A method as claimed in any of the preceding claims, 
30 further including: 

applying an additional refocusing RF pulse af- 
ter the spin echo and the first and second gra- 
dient echoes to induce a second spin echo at 
35 a time when the first and second species are 

again in-phase; and, 

generating a pair of additional gradient echoes 
after the second refocusing pulse at times when 
the first and second species are 180° out-of- 
40 phase. 

7. A method as claimed in any of the preceding claims, 
further including phase encoding the spin and gra- 
dient echoes with a common phase encoding in 

45 each of a plurality of repetitions. 

8. A magnetic resonance imaging apparatus which in- 
cludes a primary magnet (12) for generating a tem- 
porally constant magnetic field, a gradient magnetic 

50 field coil (22) for generating magnetic field gradients 
across an imaging region, one or more radio fre- 
quency coils (26) for inducing and refocusing mag- 
netic resonance of dipoles of first and second spe- 
cies within the imaging region and for receiving 

55 magnetic resonance signals from the resonating 
first and second species dipoles in the examination 
region, the first and second species dipoles having 
resonance frequencies which differ by a frequency 
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difference Scd, characterized by: 

a timing and control circuit (40, 42) for control- 
ling the radio frequency and gradient coils to (1 ) 
excite magnetic resonance, (2) refocus the 5 
magnetic resonance to induce a spin echo (54), 
and (3) induce gradient echoes (64, 68) cen- 
tered at times nn/5co before and after a center 
of the spin echo, where n is an odd integer; 
a reconstruction processor (86) for (1) recon- 10 
structing magnetic resonance signals from the 
spin echo into a spin echo image, (2) recon- 
structing magnetic resonance signals from the 
first gradient echo into a first gradient echo im- 
age, and (3) reconstructing magnetic reso- 1$ 
nance signals from the second gradient echo 
into a second gradient echo image. 

A magnetic resonance imaging apparatus as 
claimed in claim 8, further including: 20 

a phase map generator (90) which generates a 
phase correction map from the spin echo image 
and the first and second gradient echo images; 
a circuit (110, 116) which corrects one of the 25 
gradient echo images with the phase map to 
generate a phase map corrected gradient echo 
image; 

an image adder (118) which adds the spin echo 
and the phase corrected gradient echo images 30 
to generate a first species dipole image (112); 
an image subtractor (120) which subtractively 
combines the spin echo and the phase correct- 
ed gradient echo images to generate a second 
species dipole image (114). 35 



55 




7 



EP 0 745 865 A1 




8 



EP 0 745 865 A1 




9 



EP 0 745 865 A1 



European P^ot EUROPEAN SEARCH REPORT ***** Nv ** 

office EP 96 30 33G0 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation of docomo* whti indication, where appropriate, 
of relevant passages 



Relevant 

tu claim 



CLASSIFICATION OF THE 
APPLICATION (InLCL6) 



D»A 



MAGNETIC RESONANCE IN MEDICINE, 

vol. 18, 1991, 

pages 371-383, XP0G82O9847 

G.H. GLOVER, E. SCHNEIDER: "Three-Point 

Dixon Technique for True Water/Fat 

Decomposition with B0 Inhomogenei ty 

Correction 0 

* pages 371 - 377: paragraphs I and II * 

* pages 379 and 380: paragraph IV * 

US-A-5 051 699 (M. HANAWA, K. NAKABAYASH I ) 

* column 1, line 6 - line 20 * 

* column 2, line 61 - column 3, line 24 * 

* column 4, line 26 - line 68 * 

* column 6, line 40 - column 7, line 39 * 

* figures 3>6»7 * 

US-A-5 250 899 (J. L1STERUD, T. CHAN) 

* column 3, line 21 - column 4, line 52; 
figure 1 * 

EP-A-0 417 284 (YOKOGAWA MEDICAL SYSTEMS, 
LTD) 

* page 4, line 46 - page 6, line 29; 
figure 1 * 



The present search report has been drawn up far ail claims 



1-9 



G01R33/485 



1-3,5,7 
8 



1.8 



1-3,5, 
7-9 



TECHNICAL FIELDS 
SEARCHED (Inta.6) 



G01R 



Plac* af itarcb 



THE HAGUE 



11 September 1996 



Volmer, W 



CATEGORY OF CITED DOCUMENTS 

X : particularly relevant if taken alone 

Y : particularly relevant If combined with another 

document of the same category 
A : technological background 
O : noa-writtca disclosure 
P : Intermediate document 



T : theory or principle underlying the invention 
E : earlier patent document, but published on, or 

after the Tiling date 
D : document dud In the application 
L : document cited for other reasons 

& : member of the same patent family, corresponding 
document 



10 



